Integrability of two-component nonautonomous nonlinear Schrodinger equation 
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We investigate the integrability of generalized nonautonomous nonlinear Schrodinger (NLS) 
equations governing the dynamics of the single- and double-component Bose-Einstein condensates 
(BECs). The integrability conditions obtained indicate that the existence of the nonautonomous 
soliton is due to the balance between the different competition features: the kinetic energy (disper- 
sion) versus the harmonic external potential applied and the dispersion versus the nonlinearity. In 
the double-component case, it includes all possible different combinations between the dispersion 
and nonlinearity involving intra- and inter-interactions. This result shows that the nonautonomous 
soliton has the same physical origin as the canonical one, which clarifies the nature of the nonau- 
tonomous soliton. Finally, we also discuss the dynamics of two-component BEC by controlling the 
relevant experimental parameters. 
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I. INTRODUCTION 

The soliton is ubiquitous in nature and its existence 
is due to the balance of two physical features, the dis- 
persion and the nonlinearity in a nonlinear system, as 
observed by Zabusky and Kruskal The particle-like 
nature of the soliton has attracted much attention in the 
last decades since its potential application in optical soli- 
ton communication 0, H[ and its rich dynamics shown 
in different systems. Particularly, the realization of the 
Bose-Einstein condensate (BEC) in 1995 d, [f| provided 
an ideal ground to study the dynamics of the soliton 
by using the Feshbach resonance 0, 0], a technique to 
change the scattering length of atoms, and by tuning the 
harmonic external potential applied. This system can 
be well described by the so-called nonautonomous non- 
linear Schrodinger (NLS) equation with the harmonic ex- 
ternal potential. The dynamics of the soliton governed by 
this equation shows quite complicated behavior in com- 
parison with its canonical counterpart [![. While some 
early studies |8l-tl0l| have already broken the concept of 
the canonical soliton, the study of exact soliton-like so- 
lutions of nonautonomous NLS equation and discussion 
of their nature have attracted much attention in recent 
years [TM^. 

After the systematic ways to find the exact soliton-like 
solutions of the nonautonomous NLS equation have been 
proposed [HI [l2| , the attention has been focused on the 
nature of the solutions obtained. In 2003, Kruglov et 
al. argued that the self-similar solutions they obtained 
were not soliton solutions because of the non-integrability 
of the equation. However, under certain conditions the 
nonautonomous NLS equation is integrable [HI, UH, [l7l - 
[l9j . Nevertheless, these soliton-like solutions are quite 
different from the canonical ones, they move with vary- 
ing amplitudes and speeds. Thus Serkin, Hasegawa and 



Belyaeva [15| proposed the concept of nonautonomous 
soliton. Thus a question arises naturally: Could the 
nonautonomous and canonical solitons be viewed in a 
unified way? 

In this paper we indicate that the unified view is actu- 
ally the concept of balance between different features. 
First of all, we consider the single-component Gross- 
Pitaevskii (GP) equation with the time-dependent dis- 
persion, nonlinearity and harmonic external potential. 
Although the integrability condition of this equation has 
been obtained in the literature (IHIlT^ . its physical mean- 
ings has not been discussed in a deep way. Here we point 
out that in the presence of harmonic external potential 
such as BEC, the stability of the nonautonomous soliton 
can be attributed to two kinds of competitions and their 
balance. One is the competition between the nonlinearity 
and the dispersion, which is the standard one, and the 
other is the competition between the kinetic energy (dis- 
persion) and the harmonic external potential. The latter 
one has not been pointed out explicitly in the literature. 
To further deepen understanding of the picture, we con- 
sider two-component GP equation, which includes the 
intra- and inter-interactions. The procedure obtaining 
the integrability condition by using the Painleve analysis 
[2p| in this case is proved to be nontrivial but the result 
is physically straightforward, namely, for each compo- 
nent the competition between the kinetic and potential 
energies balance the competitions between all possible 
combinations of the dispersion and nonlinearity includ- 
ing intra- and inter-components (see below). This re- 
sult shows sufficiently that nonautonomous soliton and 
its canonical counterpart can be viewed in a unified way 
from the point of view of balance. Finally, based on the 
integrability condition in two-component case we explore 
the dynamics of nonautonomous two solitons and the spa- 
tial separation of binary BECs in both integrable and 
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non-integrable cases. 

This paper is organized as follows. In Section II we re- 
formulate the integrability condition of single-component 
nonautonomous NLS case. In Section III we explore the 
integrability for the double-component case. The inte- 
grability condition is composed of the combinations of 
different competition features, in which the basic form is 
similar to the single-component case. In Section IV we 
consider some applications in obtaining nonautonomous 
soliton solution of the two-component nonautonomous 
NLS equation and phase separation of binary BECs. Fi- 
nally, Section V is devoted to a brief summary. 



II. INTEGRABILITY CONDITION OF 
SINGLE-COMPONENT CASE 



Below we first consider the one-dimensional nonau- 
tonomous NLS equation 

d d 2 
i—u(x,t) + ef(t)—u(x,t) + 5g(t)\u(x,t)\ 2 u(x,t) 

~V(t)x 2 u(x,t) = 0, (1) 

where e and S are dispersion and nonlinearity coefficients 
in the standard NLS equation and f(t) and g(t) are 
dispersion and nonlinearity managements, respectively. 
V(t) is the time-dependent harmonic external potential. 
Eq. ([I]) has been extensively studied in the literature by 
various methods like the Lax pair (Til [Til [2l| and the 
similarity transformation [T^]. In Ref. |l7| . we used the 
Painleve analysis to explore the integrability condition of 
Eq. (TTJ) and proposed a systematic way to find the exact 
solutions of Eq. (TJJ directly from the known solutions 
of the standard NLS equation. These solutions can be 
obtained only under the integrability condition 

Eq. ((2) is proved to be equivalent to that in Refs. [lj^ 
figj l . In the absence of the harmonic external potential, 
Eq. (J2J is reduced to the balance condition between the 
dispersion and nonlinearity. In this case, one obtains 



g(t) 



/(*) 



ci+c 2 J f(t')dt' : 



(3) 



where c\ and c 2 are two constants determined by the ini- 
tial dispersion and nonlinearity. If one takes g(0) = 1 
and g'(0) = 0, then a = /(0) and c 2 = /'(0)//(0). Eq. 
© clearly manifests the balance between the dispersion 
and nonlinearity. In contrast to the canonical case both 
the dispersion and nonlinearity can be time-dependent 
but the system is still integrable and the soliton keeps 
stable while its shape and velocity change, as shown in 
Ref. [r|[T7]]. When the external potential is present, as 
in the case of BEC, it can be balanced by the kinetic en- 
ergy (dispersion term) , as shown in the first term of Eq. 



@. As a result, Eq. © can be viewed as the competi- 
tions between the kinetic energy and external potential 
and between the dispersion and nonlinearity. This can 
be further demonstrated by considering two-component 
nonautonomous NLS equation below. 



III. INTEGRABILITY CONDITION OF 
TWO-COMPONENT CASE 



The balance can be further clarified by considering 
two-component NLS equation 



d d 2 
i—Ui(x,t) + Eifi(t)-g-^Ui(x,t) + Sugu(t)\ui(x,t)\ Ui(x,t) 

+S i:j gi : j(t)\u : j(x,t)\ 2 u l (x,t) - -Vi(t)x 2 Ui{x,t) = 0, (4) 



where the subscripts i = 1, 2 and j = 3 — i denote the two 
components, respectively. As in the single component 
case, fi(t) and gu(t) represent the time-dependent disper- 
sion and nonlinearity in the intra-component and Vi(t) 
is the harmonic external potential. Besides the intra- 
component interaction, there is also inter-component in- 
teraction gij{t). 

To explore the integrability of Eq. ([4]) we follow the 
Painleve analysis [2(|. First of all one can expand Ui(x, t) 
and their complex conjugate Vi(x,t) = u*(x,t) on a non- 
characteristic singularity manifold ip(x,t) = x — <p(t) as 
follows 



m(x,t) = ^2u l , k (t)ip(x,t) k \ 

k=0 

oo 

Vi(x,t) =^2v itk (t)^(x,t) k - 1 . 



(5) 
(6) 



fe=0 



Inserting Eqs. ([S]) and ([5]) into Eq. (jU), performing the 
standard leading-order analysis and studying the polyno- 
mial of (f>(t), one obtain the compatibility condition for 
the two-component NLS equation as follows 



X{t)+{2a 11 {t)+2a 22 {t)-a 12 (t)-a 21 (t)) \ -lng(t)\ = 0, 

(7) 

wher ^ 9® = frfSfeg' ou(t) - *g=? and a 12 {t) = 
— T77Y ■ a 22 (t) and 021 it) can be obtained by exchanging 
the indices of 1 and 2 in the expressions of an(t) and 
o-i 2 (t). x(t) is complicated but there is a simple picture, 
it consists of the balances between the kinetic energy and 
potential energy and the possible pairs of dispersions and 
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(a): standard one 



(b): single-comp. 



f**V 



+ 



(c): two-comp. 
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FIG. 1: The unified physical picture for the balance between 
the different competitions for (a) the conventional NLS equa- 
tion, (b) the single-component NLS equation and (c) the two- 
component NLS equation i = 1, 2 (j = 3 — i) with harmonic 
external potential. 



respectively. By the similar way, the £ functions of the 
component 2 can be obtained by inter-exchanging the 
indices of 1 and 2 in the above four functions. These 
equations uncover the basic features for the integrabil- 
ity condition of the two-component NLS equation, as in 
the single-component case. That is the balance between 
the competition between the kinetic energy and the har- 
monic external potential and the competition between 
the dispersion and nonlinearity. In the two-component 
case, for each component there is one pair for the ki- 
netic energy and the harmonic external potential, which 
competes with the different pairs of the dispersion and 
the nonlinearity, namely, f\ f> jn, a contribution from 
the intra-component; /j f-> 512, a contribution from the 
inter-component; f 2 <-> 321 and f 2 <-> 322, due to the 
existence of the component 2 and relevant competitions 
between the dispersion and the nonlinearity in it. To clar- 
ify the above features, in Fig. [1] we show schematically 
this picture and make comparison with the cases in the 
single-component and the standard NLS equations. This 
result indicates that the integrability condition in these 
equations can be viewed as a unified picture, irrespective 
of the standard or nonautonomous systems. 



nonlinearities in the two-component system, namely, 
X(t) = a 11 (t)' ' 



lUl,/ 1 n 9ll (t) + (1^2) 



+a u (t)(Xg(t) - 1)6,Ab Ji3 W + (1 ^ 2) 



+ aia(*) - 



dl 2 (t) 



a 12 (t)Xg(t) 



£i,/ 3 o ff21 (*) + (1^2) 



£i,/ 2 o 9M (*) + (1^2).(8) 



Here the expression of each (1 f± 2) is obtained by inter- 
changing the subscripts of 1 and 2 of the left term given in 
each row correspondingly in the right of the above equa- 
tion, and A = 3773^7- Four £ functions for the component 



1 read 



(t) = 2e 1 f 1 mit) + ^ In ' ' 



dt 2 gil (t) 



d, m d h{t) 



(9) 



Zi, h „ g „(t) = 2e l f l (t)V 1 (t) + 4ln g22 i!| h{ % 

dt 2 021 (*) 912 (*) 



d g 2 i(t) d g 22 (t) fi{t) 

— In ttt 9\2(t)— In- 



dt g 22 {t)' dt g 2 i{t) g 12 (ty 

£i, A++flB1 (t) = 2e 1 f 1 (t)V 1 (t) + ^In^- 

d fi{t) d f 2 (t) 

-—In -r-7-r52l (*) -77 In ■ 



(10) 



6,/ 2 « S22 (t)-2 £l / 1 (^(i) + ^ln ffl2f/i 



(11) 



dt 2 9ii(t)g 22 (t) 

.^mi^^^i^m (12) 

<a / 2 (0.9i 2 W * 511 (*) 522 (t) k ; 



IV. APPLICATIONS 

The generalized integrability condition of Eq. ([7]) for 
the two-component nonautonomous NLS equation has 
not been reported in the literature. It can be used to 
investigate unique phenomenons in many practical two- 
component nonlinear systems. To further explore it, first 
we consider the symmetric case, such as e± = e 2 = e, 
fi(t) = h{t), gu(t) = 922 (t) = gv»(t) = ff 2i (t) = g(t) 
and Vi(t) = V 2 (t) = V(t). Thus Eq. © reduces to 

^wu+i^-i^u^w)' - (13) 

which is the same as that in the single-component case. 
Furthermore, if e = 1/2 and fit) = 1, Eq. (ITBl gives a 
constraint condition on the applied external potential 



V(t) 



1 d 2 
g(t) dt 2 



g(t) 



g(t) 2 \dt 



-At) 



(14) 



This expression is completely consistent with those re- 
ported in the literature 22] by using different methods. 
Under this condition, one may take a general transforma- 
tion to obtain exact solution of the coupled NLS equation 
flU), as in the single-component case [lj| 

Ui (x,t) = U i (X(x,t),T(t))e i < x ^- i < t \ (15) 

where Ui(i = 1,2) satisfy the standard coupled NLS 
equation 

+8 l] \U ] (X,T)\ 2 U l {X,T)= 0.(16) 
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The explicit forms of X(x,t),T(t),a(x,i) and c(t) 
can be obtained analytically as follows a(x, t) — 

Ii7(t) {it ln 1$) x2 + C W) X ~ eG2 1 7W dt ' X ( x ' ^ = 
^-^Cf^dt, c(t) = |ln(f|), and T(t) = 

/ ym~dt, where C is an arbitrary constant. Eq. ([To]) has 

a standard two-soliton solutions [23[ , from which the cor- 
responding nonautonomous two-soliton solutions of Eq. 
((4]) can be explicitly obtained by Eq. (TTSl) . To understand 
the influence of Feshbach resonance on nonautonomous 
two-dark-bright soliton behavior under the integrability 
condition Eq. (TlUl) , we investigate the dynamics of these 
solitons with Feshbach management. These solitons can 
be realized in a two-component system such as different 
hyperfine spin states [24| . g(t) is modulated periodically 
as g(t) = 1 + Asm(a;t), where A and lj are the relevant 
modulation parameters. The corresponding trap poten- 
tial is tuned by Eq. (jT4"]) . In Fig. [5] we show the nonau- 
tonomous two-dark-bright solitons dynamics obtained (c, 
d) and also compare them with the canonical ones (a, b) . 
The result shows clearly that these solitons remain sta- 
ble, as the conventional ones, but their shapes and speeds 
change with time. It is significant to guide experimental 
control over the dynamic behavior of the muti-solitons in 
binary BECs. 
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FIG. 3: (Color online) Spatial separation of two components 
in the nonintegrable (a, c) case for A = 1 and the integrable 
(b, d) case for A = 0.5 and uj = 2. 



Unfortunately, the explicit solution of Eq. Q can not be 
obtained directly from that of Eq. (fTrJ|) in this case, but 
it is enough to use the Gaussian wave packets to explore 
the spatial separation problem. They read [251 ] 



a 




FIG. 2: (Color online) The standard two dark-bright solitons 
(a, b) and the corresponding nonautonomous dark-bright soli- 
tons with the Feshbach management (c, d). The parameters 
used are A = 0.5 and cj = 1. The other parameters used are 
oi = 1, a 2 = 0.77, bi = b 2 = 0, c = and r = 1. 

Finally, we consider the phase separation phenomenon 
in binary BEC. For simplicity, we let Sn = 822 = <5i2 = 
fei = -1, hit) = 1, £ S : = 1/2,' 9u(t) = 1 and V t {t) = V{t) 
(i = 1, 2). Thus Eq. © becomes 

(17) 



U X {x,t) = A(t)e~ {j ^$- ■ e i{C{t)+D(t)x+E(t)x*) ( 
U 2 ( X ,t) = A(t)e~ ( »W e i{C(t)-D(t)x+E(t)x 2 ) ^ 

where A(t), B(t), C(t), D(t), E(t) and W(t) are the vari- 
ational parameters and are determined b y a set of ordi- 
nary differential equations given by Ref . [25| ■ 

We solve this set of equations in two different cases. 
The first one is to increase intercomponent interaction 
linearly but fix the harmonic external potential as shown 
in Fig. [3Ja). Initially the two components are mixed and 
keep miscible for some time. Subsequently the spatial 
separation happens. This is in agreement with the pre- 
vious observation [25|. However, it is necessary to point 
out that this case is not integrable. In contrast, in the 
second case we consider the integrable constraint on the 
harmonic external potential given by Eq. (|17p . Here the 
intercomponent interaction assumed varies around 1 pe- 
riodically, as shown in Fig. EJb). In this case the applied 
external potential is modulated periodically following the 
intercomponent interaction. The two components in this 
case are found immiscible, even if the miscible condition 
<7i2 (t) < 1 is satisfied frequently. This is a consequence of 
the change in the external potential, or rather the pres- 
ence of the potential barrier (V(t) < 0) in some time, 
which leads to the two-component separation shown in 
Fig. ED^d). The two different results provide possible ways 
to realize two-component spatial separations in binary 
BECs. 
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V. SUMMARY 

By analyzing the integrability conditions of the nonau- 
tonomous NLS equations we found that the nature of the 
nonautonomous soliton is due to the balance between 
different competition features: the kinetic energy versus 
potential energy and the dispersion versus nonlincarity. 
This picture is consistent with the conventional one which 
originates from the balance between the dispersion and 
nonlincarity. On the other hand, the integrability condi- 
tions in two-component case is useful to guide experimen- 



tal control over the dynamic of nonautonomous solitons 
and two-component separation in binary BECs. 
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